A variational approach for calculating Auger electron spectra: going beyond the 
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We propose a novel variational method to calculate the two-hole propagators relevant for Auger 
spectroscopy, which can be thought of as an intermediary step between the full solution (which is 
difficult to generalize to systems with partially filled bands) and the impurity approximation. Like 
the former, our solution has the full translational invariance, and like the latter, it can be generalized 
to systems with partially filled bands. Here, we compare both our variational approximation and the 
impurity approximation against the exact solution for a simple one-dimensional model with filled 
bands, and show that when both are valid, the variational approximation is superior. 



I. INTRODUCTION 

Spectroscopic measurements are a powerful set of tools 
for probing various aspects of many-body physicsJ^ 
Among these, Auger electron spectroscopy (AES) pro- 
vides information about the local atomic multiplet 
structure, on-site interaction strengths and the crystal 
fieldsi^"— In transition metal oxides, the AES of the tran- 
sition element can be supplemented by the O KLL Auger 
spectra, resulting in additional information about the O 
on-site repulsion energy as well as interactions between 
holes located on nearest-neighbor transition metal and O 
ionsi^ Such information is vital for understanding corre- 
lated materials, which is why AES has been the subject 
of sustained research for a long time. 

The Auger process consists of the decay of a core hole 
into two final state holes plus an Auger electron, and is 
mediated by onsite Coulomb interactions. The two fi- 
nal state holes are initially located on the atomic site 
where the original core hole was created by the high en- 
ergy X-ray. One of the most studied cases has these 
two final state holes residing in the valence band, and 
goes by the name of Core- Valence- Valence (CVV) Auger 
spectroscopy. From a theoretical point of view, the eas- 
iest case to handle has a full valence band except for 
the two Auger holes. The central quantity of interest 
is the two hole Green's function, which for an otherwise 
full band can be calculated within the two-step approxi- 
mation using the Cini-Sawatzky theoryii"— The resulting 
two hole spectral function, multiplied by momentum de- 
pendent matrix element factors, provide the theoretical 
predictions for AESJ^ Many extensions of this approach 
have been carried out to incorporate various aspects such 
as dynamical screeningfiir— off-site interactions^^ over- 
lap effectsr^"— and one step formulationii^ We refer the 
reader to a recent reviewiS for more details. These ef- 
forts have led to spectacular success in explaining AES 
for materials such as Cu and CU2O. 

However, the problem of understanding AES for sys- 
tems with a partially filled valence band, for example 
the oxides of transition elements Ni, Co, Fe, Mn, etc.. 



remains open, because the two hole spectral function is 
very challenging to compute in this case. Limited suc- 
cess has been achieved employing variants of the bare 
ladder approximations^^ and assuming low hole density 
in the bands. A completely different approach is to use 
the Anderson impurity approximation, whose underlying 
idea is as follows: in transition metal oxides, the transi- 
tion metal atoms are typically connected to each other 
through oxygen ligands. The simplest example in one 
dimension is sketched in Fig. [Ija), and has transition 
metal atoms intercalated with O atoms. In the impurity 



approximatiouf 



I the full problem is simplified to that 



of a single transition metal atom coupled to the bath of O, 
as sketched in Fig.[ljb). This greatly simplifies the calcu- 
lation and is a reasonable step towards understanding lo- 
cal multiplet structures. However, because the symmetry 
of the problem is lowered, momentum-resolved spectral 
weights cannot be calculated. 

Here, we propose a variational approach for finding the 
two-hole Green's functions needed for AES, which can be 
thought of as an intermediary step between the impurity 
approximation and the full lattice calculation. Like the 
former, it has a reduced Hilbert space and can therefore 
be generalized quite easily to systems with partially filled 
bands. Like the latter, it has full translational invariance 
so that momentum-resolved spectral weights can be cal- 
culated. We argue that the location of the spectral weight 
components ignored within our approximation can be in- 
ferred apriorily but should be small for oxides and other 
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FIG. 1. (a) Sketch of the ID periodic Anderson model. 'TM' 
and 'O' represent transition metal and O ions, respectively, 
(b) In the impurity approximation, an impurity 'TM' ion is 
coupled to the ID 'bath' of O. 
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materials where AES is a useful probe. Moreover, the 
variational space can be systematically enlarged to ver- 
ify the relevance of some of the neglected terms. 

In this paper we present the underlying idea and the 
general formalism of this variational approach, which is 
based on a recently developed method to calculate many- 
particle Green's functions)^"— Here, we apply it to the 
simple model sketched in Fig. [T] and assume otherwise 
full bands, so that we can benchmark it against the ex- 
act solution available in this case. For completness, we 
also show impurity approximation results for this model. 
This allows us to gauge the advantages and disadvan- 
tages of both these approximations and to understand 
for what regions in the parameter space our variational 
approximation is valid. This information will guide us, 
in future work, to using it for appropriate systems where 
the transition elements have partially filled d orbitals. 

The paper is organized as follows. In Sections II we 
define the Hamiltonians and in Section III we discuss the 
methods to calculate two-hole propagators exactly and 
with the two approximations. In Section IV we present 
the exact results and compare them with those predicted 
by our variational approximation as well as the impurity 
approximation for two different topologies of coupling be- 
tween the impurity and the bath of O. Section V contains 
the summary and conclusions. Some of the details are 
presented in the two Appendixes. 

II. HAMILTONIAN 

A. The periodic Anderson model 

The periodic ID Anderson model shown in Fig. [IJa), 
is defined by the Hamiltonian: 

"H-p = Htm + Ho + Hhyb 

where Htm describes the chain of TM atoms, which 
for simplicity are assumed to have only two (spin- 
degenerate) states each, and is given by: 

i i,{7 

The O "bath" is described by a ID Hubbard model: 

Ho = -t ^{cl^^Ci+i^a + h.c.) + Uo'^ ni,t'^^,l 

i,a i 

while the TM-0 hybridization is described by 

Hhyb = ^^X! ['^l'^ + + ■ 

i.cr 

Here, all creation operators are electron creation oper- 
ators, with CiCT and di^^r the operators for O and TM 
orbitals, respectively, with the convention that i*'* TM 
atom is placed to the left of the i"* O atom. As usual. 



nd,i,a = d]^dia,'riia = cJ^Cicr. For book-keeping purposes 
we assume that there are TV TM and N O sites, but in 
the final calculation we let A^ oo. We also assume that 
only the 2p O orbitals pointing towards the TM neigh- 
bors are relevant to the calculation. Generalization to 
models that include more O and/or TM orbitals, as well 
as lattices in higher dimensions, arc discussed below. 

For AES in a system with full bands, we start from 
the completely full ground state \il)p = Hi cr '^I a*^! cr|0) 
of energy = N{2A + U + Uq), and create two 
holes. Because the periodic Hamiltonian is invariant 
to translations and because we want to describe states 
where the two holes could be at the same TM site, 
we choose two-hole basis states with a total momen- 
tum k and zero spin projection^ namely \k,n,dd) ~ 
_^^^gifc(fl,-Hna/2)^^ ,^^^^^^|j^^p if both holes arc on 

TM sites, \k,n,dc) = ^ e'''''^"+""^'''^^tc^+n,i|f^)p 
and \k,n,cd) = e*'=(^'+"'^/2)c,,^d,+„,;|f})p if one 

hole is on a TM site and the other is at an O site, 
and finally \k,n,cc) = e'Mfl.+W2)c,_^c,+„,;|f})p 
if both holes are in the O bath. Here a is the lattice 
constant and n = — ^ + !,■•■, takes all possible val- 
ues consistent with the periodic boundary conditions. 
Taken together, these states constitute a full basis for the 
Hilbert subspace containing states with total momentum 
k and zero spin projection. 

The aim is to find the propagator 

Gddik,0,u) = {k,0,dd\Gp{u)\k,0,dd) 

where Gp{uj) = [uj + ir] — {Hv — ^n)]~"^ the re- 
solvent for Hp, because its associated spectral weight 
Add{k, 0,io) = — ilm[Gdd(fc, 0, cj)], which has poles at en- 
ergies UJ = E2h{k) — for any two-hole eigenstate with 
total momentum k and energy E2h{k), is proportional 
to the momentum-resolved spectral intensity measured 
by AES. As discussed below, our solution also provides 
the values of many other propagators beside Gdd{k, 0, uj), 
from which other useful information can be gleaned. 

B. Anderson impurity problem 

We reduce the full periodic problem to the Anderson 
impurity problem sketched in Fig. [Ijb), by simplifying 
the Hamiltonian to: 

Hi = Htm,i + Ho + Hhyb,! 

where 

Htm, I = Und,^nd,i + A ^ Ud^a 

a 

describes the impurity TM site, the O bath is described 
by Ho as before, and 

^^iiyb,i = -V^^{d\cx.a + rf^-cco- + h-c\ 



3 



is their hybridization, with the impurity taken to be lo- 
cated between the bath O labeled '0' and '1'. Here, again, 
all operators are electron operators and dj. is the creation 
operator for the orbital of the TM impurity site. 

The filled-band ground state is now = 

Tier '^o- Hi cr|0) ^iid its Corresponding energy is = 
2A+U+NU0. To calculate AES-relevant spectra we con- 
sider two-hole excitations in this ground state by remov- 
ing two electrons with opposite spins. Since invariance to 
translations is lost, the generic real space states of inter- 
est arc now \icr,ja') = Cj^c^-^, \d(j,ia') = d^c^^,\n)j 
and \dd) — d^di\n)j. Now we are primarily interested 
in calculating the two-hole impurity Green's function 

I 



Gdd{^) = {dd\G I {uj)\dd) where Gi{uj) = [ui + ir]— ('Hi — 
E[^)]^^, and its corresponding two- hole spectral function 
Addiuj) -ilm[Gdd(w)]. 



III. METHOD 

Wc begin with the exact solution for the periodic An- 
derson model. To find the propagator of interest to us, 
Gdd{k,0,u!), we generate its equation of motion (EOM) 
from the identity: Gp{uj){uj — 'Hp + Ef^ + irj)=I. Calcu- 
lating its diagonal matrix element for \k,0,dd), we find: 



{UJ + 2A + U + i'n)Gdd{k, 0,Lu) = l + V e^Gcdik, 1, uj) + Gcd{k, 0, uj) + e^Gddk, -1,lu) + Gddk, 0,lu) 



where Gap{k, n, uj) = {k, 0, dd\Gp{u!)\k, n, aj3) for a,j3 = 
c, d. In other words, because the Hamiltonian links the 
state \k,0,dd) to states with one hole on a TM site and 
one on a neighboring O, the EOM links Gdd{k,0,uj) to 
propagators corresponding to such states. Their EOM 
can be generated similarly, and we obtain an infinite se- 
quence of coupled linear equations. 

To solve it, we couple the propagators with holes at 
the same distance n in a vector: 



14 = 



(Gdd{k,n,uj) 
Gdc{k,n,uj) 
Gcd{k,n,uj) 
Gcc{k,n,uj) 



and note that for any n 7^ 0, the EOM can be grouped 
in the simple recurrence relation: 

for any given k and uj. Here 7^, /3n and a„ are simple 
4x4 matrices that are read off directly from the EOM. 
We note that one can always group the EOM in such 
simple recurrence relations, even for models which allow 
longer-range hopping^^ and/or in higher dimcnsionsj^ 
For n = 0, the recurrence relation also has an inhomoge- 
neous term: 

joVq = X + /3oVi + aoV-^i, 

where X'^ — (1, 0, 0, 0) for this problem. The solution of 
such recurrence relations has been discussed extensively 
elsewhere 3^"— Briefly, we must have V^i -t- as |rt| — >■ cxj, 
because the Fourier transform of these propagators are 
the amplitudes of probability to have the two holes evolve 
from being on the same TM site to being n sites away 
from each other, in a given time. As \n\ 00 this be- 
comes very unlikely, and the presence of the broadening 
77 which introduces an artificial lifetime 1 /rj makes it even 



I 

less so. As a result, for rt > 1 we have Vn = An{k, uj)Vn-i 
where An = [7,1 — f5nAn+i]~^an is calculated starting 
with Am = for a sufficiently large cutoff M. Sim- 
ilarly, for n < — 1 wc have Vn = Bn{k,LL))Vn+i where 
Bn = [7n — UnAn-i]^^ Pn is Calculated starting with 
B_\M\ = 0. Using Vi = AiVq and V^i B^iVq in 
the n = equation gives: 

= [70 - PqAi {k,uj)- aoB^iik, ij)]-^ X. 

This gives us Gdd{k,0,uj) as the top entry in Vq. All 
other n = propagators, as well as those with |7i| < Af , 
can also be then calculated efficiently. Projecting on a 
different state than (fc, 0, dd\ also simply requires using a 
different X, so other propagators can also be found easily 
with this scheme. 

In principle, this method generalizes straightforwardly 
to lattices in any dimension and with any topology, so 
no approximations should be necessary. In practice, 
however, the computational cost quickly becomes pro- 
hibitive. For instance, for models with nearest-neighbor 
hopping in higher dimensions, one must group together 
in Vn all propagators where the holes are separated by 
nxaex+nyaey + . . . with |na;|-|-|ny|-|-- •• = rt, i.e. with the 
same Manhattan distance^i As a result, the dimension 
of Vn increase roughly like 7i'*~^. Adding more orbitals 
at the TM/0 sites will further amplify the problem, by 
increasing in a combinatorial fashion the number of pos- 
sible propagators with the same {ux, Uy, . . .) separation. 
As a result, while the solution can still be cast in terms of 
continued fractions of matrices, their dimensions increase 
fast with n resulting in significant computational costs. 
This is why efficient approximations are needed. 

As already mentioned, one much-employed option is 
the impurity approximation. Its solution for Gdd{^) can 
also be formulated in terms of continued fractions of ma- 
trices (we present the details in the Appendix, since we 
are not aware of a prior similar solution for this prob- 
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lem) . This approximation reduces the number of possible 
propagators by removing aU but one TM site from the 
problem. As a result, there is now only one state with 
both holes at the impurity TM site as opposed to N"^ in 
the full periodic problem (i.e., including all allowed total 
momenta), and only 2N combinations with one hole at 
the TM and one at an O site, as opposed to 2iV^ options 
in the full problem. The number of states associated with 
both holes in the O bath is not changed. 

The total number of distinct propagators is therefore 
reduced from 47V^ in the full problem to {N + 1)^ in the 
impurity limit, suggesting that the latter is considerably 
more efficient. However, the full problem can be solved 
individually for each of the allowed N total momenta, 
since the translational invariance guarantees that prop- 
agators with different momenta do not mix. Thus, one 
needs to solve N distinct problems with AN propagators 
each, as discussed above. In the impurity problem, loss 
of translational invariance means that all propagators are 
coupled to one another through the EOM. From this per- 
spective, it is far less clear that using the impurity ap- 
proximation is computationally more efficient, although 
to fully settle this one needs to take into account of all 
other possible symmetries and the effects of the trunca- 
tion (for large system, the cutoff M <C N/2). In practice 
the O bath is often replaced by a featureless density of 
states (we do not make this further approximation in our 
calculations). This certainly makes the impurity problem 
very efficient, but it also looses all information regarding 
the role of the O bath topology, on top of the loss of 
ability to momentum-resolve the AES spectral weights. 

From this analysis, it is clear that a better strategy 
would be to lower the total number of propagators while 
maintaining translational invariance. This is the basis 
for our proposed variational method. Our proposal is to 
remove the propagators which have the holes at different 
TM sites. This is equivalent with excluding these states 
from the Hilbert space, which is why this is a variational 
approximation. Mathematically, this is easily achieved 
by setting Gdd{k, n, w) = for all n 7^ in the EOM dis- 
cussed above, leading to vectors y„, n 7^ of dimension 
3 instead of 4. Of course, the saving is not big for this 
simple problem, but it starts to become considerable in 
higher dimensions and/or for more TM/0 orbitals. 

Physically, what this means is that both holes start at 
the same (any) TM site. Eventually both may hop into 
the O bath and move through the system, but when- 
ever they happen to hop back to TM sites they must 
both go to the same TM ion. For systems with strong 



correlations, i.e. where the states with the holes at the 
same TM site are at quite different energies from states 
with the holes as different TM sites, this should be a rea- 
sonable approach and is close to the intuitive picture of 
what happens in AES. In a way, one could view this as a 
"lattice of impurities" . Besides maintaining translational 
invariance, this variational approximation has the added 
benefit that we can infer apriorily the effect of removing 
these states, as we discuss in the next section. More- 
over, if need be one can also systematically add some 
of these states back into the calculation - for example, 
starting with the states where the holes belong to neigh- 
boring TM sites only. Their importance can therefore be 
assessed quantitatively. Generalization to systems with 
partially filled TM orbitals is also much less computa- 
tionally costly than for the full solution. More discussion 
on the meaning and relevance of this approximation is 
presented below, as we discuss its predictions. 



IV. RESULTS 
A. The periodic Anderson model 

We start by discussing the full periodic model. Be- 
sides providing the test case against which we compare 
the variational and the impurity approximations, this 
also allows us to understand the various features of these 
spectra and their dependence on the various parameters. 
Since the phenomenology which leads to the appearance 
of these various features is general, qualitatively similar 
spectra should be expected in more complex systems. 

To understand the dependence of the AES spec- 
tral weight with both holes at the same the TM site, 
^ddik,0,Lo), on the various parameters, we start by set- 
ting U = Uo = 0, t = V = 1, and varying A. 
Momentum-resolved results are shown in Fig. [2] for uni- 
formly spaced values of k from and n/a. 

To make sense of these fairly complicated spectra, 
we note that these parameters correspond to a non- 
interacting system. As a result, the two-hole spectra 
must be convolutions of single-hole spectra, which are 
easy to calculate for this simple model (note that single- 
hole spectra can also be obtained experimentally from 
Angle- Resolved Photoemission). Straightforward calcu- 
lations shows that a hole of momentum k introduced in 
the state \i})p has two possible eigenenergies: 



~U-A-Uo+2t cos(fca) ± \j {U + A ^ Uq + 2t cos{ka)f + 16V^ cos^ y 



(1) 



Roughly speaking, these correspond to either having the hole preponderantly on the TM sites (i.e., a TM-like 
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FIG. 2. (color online). AES spectral weights Add{k,0,uj) 
for the full periodic model vs. cj for k £ [0,n/a]. The curves 
are shifted vertically with increasing k. In all cases t = V — 
= f/o = 0,7? = 0.05, while (a) A = 0; (b) A = 5; (c) 
A = —5. The small vertical lines indicated the expected 
locations of band edges. See text for more details. 



band) or in a 0-like band. Note that the electron eorre- 
lations energies U, Uq enter here because the single hole 
is introduced in an otherwise full band. 

The convolution of these two one-hole continua re- 
sult in three two-hole continua, covering the ranges 



Krik) e {E\l'{k -q) + - ^ < q < ^} for 

(5,7 = ±. Their band edges are shown by small verti- 
cal markers in Fig. [2] each band has a different color. 
These expected band-edges indeed agree perfectly with 
the features seen in the AES spectral weight (a broaden- 
ing 77 = 0.05 was used in the spectral weight, accounting 
for the apparent "overflow" ) . 

For A = 0, the two upper bands overlap partially. 



3- 



3 
o 

1 - 



(a) A=5, U=5, U„=0 





-20 




-15 



(O/t 



3 
o 

* 2- 
<'° 

1 - 



1 1 1 


1 

s" 

o 
II 

-a 
< 



— 




J 






III 
(b) A=5, U= 

V_i,i 


III 
5,U„=3 

1 " 


-20 


-15 -10„/j -5 5 

















































-20 



(O/t 



FIG. 3. (color online). AES spectral weights for the full 
periodic model Add{k,0,Lj) vs. oj for k G [0,7r/a]. The curves 
are shifted vertically with increasing k. In all cases t — V — 
1, A = 5, 77 = 0.05, while (a) U = 5,Uo=0; (b) U = 5,Uo = 
3. The small vertical lines indicated the expected locations of 
band edges. See text for more details. 



and the spectral weight is distributed fairly equally be- 
tween all features. For A = ±5, most of the weight is in 
the lowest/highest continuum which contains states with 
both holes in the TM-like band and thus has the highest 
overlap with the state \k,0,dd) of interest to AES. This 
band is centered at — 2A, which is the change in energy 
if two holes are removed from TM sites if there was no 
hybridization, V = 0, and in the absence of correlations 
U = Uo = 0, and is fairly narrow since the holes cannot 
hop directly between TM sites. The middle continuum, 
with one hole in the TM-like band and one in the 0-like 
band, has less weight but is still visible. It is centered 
around —A, i.e. the energy cost for removing one elec- 
tron from the TM site, and is broader because the 0-like 
band has significant bandwidth. The third continuum, 
with both holes in the 0-like band, has very little over- 
lap with the state \k,0,dd) and is therefore only visible 
when the scale is significantly expanded, as shown in the 
insets. It is centered around the origin (for Uq = 0) and is 
the broadest of the three features because holes can hop 
directly between O sites. In the limit V ^ 0, the max- 
imum bandwidth of this continuum should be 8t, when 
k ^ Q. We see that hopping to and off the TM sites, 
controlled by V , can have a significant influence on this 
value when V ^ t. 
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While one can set any values for parameters in a theo- 
retical study, physically it makes sense to focus on cases 
where the low-energy states favor having the extra holes 
at the TM sites, so that AES is maximally sensitive to the 
TM parameters. Simple arguments, listed in Appendix 
B, show that this implies A > 2t + Uo- 

As a result, we set A = 5 and investigate the role of 
the on-site correlations. Fig. El^a) shows Add{k,0,uj) for 
U = 5,Uo = while Fig. ^h) is for C/ = 5, [/q = 3. 
Comparing Fig. ^a) with Fig. EJb), we see that the 
lowest band with the two holes in the TM-like bands has 
shifted by about 2U and is now centered around — 2(A + 
U), as expected since ~{A + U) is the energy cost for 
removing an electron from a TM site (if there was no 
hybridization with the O bath). Similarly, the middle 
band is shifted by about U to around — A — 17, while the 
upper band with the two holes in the O bath is essentially 
unchanged. Apart from these expected shifts, we sec that 
most spectral weight has moved from the lowest band 
(where it was for U ^ Q) into a new discrete peak that 
has appeared U above it. This weakly dispersing peak 
describes (anti)bound states with both holes at the same 
TM site, hence the large overlap with |fc, 0, dd). Indeed, it 
is located close to — 2A — U, where it would be expected 
to appear in the absence of hybridization with the O 
bath. For our simple model, this peak represents the 
"multiplet structure" associated with the TM element. 
For a more realistic calculation, for example for CU2O, 
this peak would be replaced by the atomic-like multiplet 
structure of the Cu as seen in experiments!^ 

Correlations at the O site have similar effects, as shown 
in Fig. [3jb). Since the cost of removing an electron from 
the O ion is now lowered hy Uo, the central band shifts 
by an additional —Uo while the upper band shifts by an 
additional —2Uo- On the other hand, the lowest band 
with the two holes in the TM-like band is essentially un- 
changed, while the "multiplet" peak is shifted to slightly 
lower energies due to level repulsion with the central con- 
tinuum. Correlations at O sites also produce an (anti)- 
bound discrete state with both holes at the same O site 
which is pushed above the upper continuum. However, it 
has very little weight and can only be seen on a greatly 
expanded scale, as shown in the inset, as a small peak 
located just above the highest-energy bandedge. 

To summarize, the location of the continua in the 
momentum-resolved AES spectrum must agree with the 
convolution of the one-hole spectra, which can be ob- 
tained from ARPES. Additional discrete peaks (or strong 
sharp resonances, if these happen to fall inside another 
continuum) indicate the presence of on-site correlations, 
and allow one to find various on-site parameters from 
the corresponding multiplet structure. Of course, be- 
cause the projection is on the |fc, 0, dd) state, the spectral 
weight is large for states which predominantly have both 
holes at the same TM site. As discussed in the previ- 
ous section, it is trivial with our method to project on 
other states, for example with both holes at the same O 
site. While the spectrum is unchanged, this will shift the 




FIG. 4. (color online). Same as in Fig. O but the AES spec- 
tral weight is obtained with the variational approximation. 



spectral weight between the different features. 



B. Variational approximation 

We now present data obtained with our variational ap- 
proximation for the same sets of parameters, in Figs. |3] 
and[5l As before, the vertical markers show the expected 
location of the band-edges, based on the convolution of 
the one-hole spectra. 

Starting with the un-correlated case in Fig. 21 we see 
that this variational approximation is very poor when 
A 0, in panel (a). While the overall spectral range 
is in fair agreement with that predicted by the two-hole 
convolutions, inside this interval there is significant dis- 
agreement. In particular, the variational calculation pre- 
dicts spectral weight at energies that should be gaped. 
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FIG. 5. (color online). Same as in Fig. [3l but the AES spec- 
tral weight is obtained with the variational approximation. 



This disagreement is not surprising, since in this case we 
are projecting out states that have energies very similar 
to the states we keep, and this is not a sensible strategy. 

As shown in panels (b) and (c), the situation improves 
when A is large enough to separate the TM-like band 
from the 0-like band. In this case, the agreement for 
the two less visible bands is reasonable, even though the 
projection may slightly over/underestimate the location 
of the band-edges and it certainly produces additional 
structure within the bands, especially the central one. 
This is not surprising since one expects a fair amount of 
hybridization between states in the central continuum, 
which have one hole in the TM-like band and one in the 
0-likc band, with the states with the two holes at differ- 
ent TM sites that arc projected out. 

Finally, the continuum with the two holes in the TM- 
like bands is replaced by discrete peaks located at roughly 
the correct energy. The disappearance of this continuum 
is expected, since we projected out precisely the states 
responsible for generating it, i.e. with holes at different 
TM sites. In the absence of correlations, the state with 
both holes at the same TM site has a similar energy, 
hence the peak. Similar conclusions hold for A = —5. 

If we turn on the correlations, this "multiplet struc- 
ture" peak moves away from the continuum. Indeed, as 
shown in both panels of Fig. [5l now the low-energy con- 
tinuum with the two holes in the TM-like band is com- 
pletely absent from the variational results but the impor- 
tant feature, i.e. the high-weight "multiplet structure" 




FIG. 6. Two ways to connect the impurity to the O bath: 
(a) trivial, and (b) non-trivial topology. The latter can be 
reduced to the former if one considers the three central ions 
enclosed by the box to form an effective "impurity" . 



peak, is essentially indistinguishable from that predicted 
by the full solution. The higher energy continua show 
some additional structure within this variational calcula- 
tion, but since this structure also appears in the absence 
of correlations there is little danger of confusing it with 
real features induced by the correlations. 

To summarize, for parameters likely relevant for the 
real materials, i.e. when the TM-like and 0-like bands 
are sufficiently well separated, the variational approxi- 
mation does a very good job in capturing the multiplet 
structure due to TM on-site correlations. One of the 
continua is projected out, but it has little weight and its 
location is apriorily known from the one-hole spectrum. 
The higher energy continua are reproduced at roughly 
the correct locations but with some additional features, 
which however also appear in the absence of correlations. 
With some care, one can use this more efficient calcula- 
tion to understand most, if not all, the important physics 
contained within these spectra. 



C. Impurity approximation 

For comparison, we also present AES spectra obtained 
with the impurity approximation. Since here the trans- 
lational invariance is lost, we can only show one Addi^) 
curve for each set of parameters. This can be roughly 
thought of as a momentum- integrated spectral weight. 

We begin by considering the simpler impurity model 
sketched in Fig. ^a.), where the TM impurity is con- 
nected to two independent semi-infinite half chains. The 
reasons for studying this case first are that (i) it is sim- 
pler to understand the resulting AES spectrum; (ii) one 
can think of the actual case of interest, shown in Fig. 
inijb), as a generalization of this case, by treating the 3 
central atoms enclosed in the box as an effective impu- 
rity with a more complicated level structure, but which 
is otherwise connected to two independent semi-infinite 
half chains; (c) comparison between these two cases will 
reveal whether the particular way in which the impurity 
connects to the O bath plays a significant role. 

Fig. [7] shows results for several sets of parameters. In 
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FIG. 7. (color online). Impurity spectral function Add{i^) vs. 
oj for the trivial topology of Fig. [6ja) with t = V = 1 and (a) 
A = f/ = (7o = 0, (b) A = 5, f/ = (7o = 0. The inset shows 
the two-hole continuum on an expanded scale, (c) A = 5, ?7 = 
2, 5, 10, Uo = 0, and (d) A = -10, -5, 2, [/ = 5, (7o = 0. 



all cases t ^ V ~ 1. Panel (a) shows the case A = 
U = Uo = 0. Since there is now no difference between 
the TM and the O sites, this describes two independent 
holes on an infinite chain with nn hopping. The spectral 
weight has the expected shape for this problem, with the 
two-hole convolution continuum ranging from [— 4t,4t]. 
This is quite different from the exact result in this limit, 
shown in Fig. [2|a,. 

In all the other panels A 7^ so that the single hole 
spectrum consists of the O band and a TM-like impu- 
rity state. As a result, the 2-holes convolution consists 
of three features. The first is a discrete state which pri- 
marily has both holes at the TM site. In the absence of 
hybridization {V = 0), this feature is located at —2A — U 
and represents the "multiplet" of primary interest. The 
second feature is a one-hole continuum describing states 
with one hole occupying the TM impurity state and the 
other moving freely in the O bath. This replaces the con- 
tinuum with one hole in the TM-likc band and one in the 
0-like band, from the exact and the variational solutions. 
In the absence of hybridization, this feature extends over 
[-A-U-Uo-2t,-A-U-Uo + 2t]. Finally, the third 
feature is a two-hole continuum with both holes in the 
O bath. This corresponds to the continuum with both 
holes in the 0-like band, but here it spans the interval 
[-2Uo - it, -2Uo + it] irrespective of the value of V. 

The results in panels (b)-(d) of Fig. [7]agree with these 
expectations. Like in the full calculation, the continuum 
with the holes in the O bath has very little weight and 
can only be seen on an expanded scale, if A 7^ 0. The 
one-hole continuum is at the expected location, while the 
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FIG. 8. (color online). Impurity spectral function Add{u}) vs. 
(jj for the non-trivial topology of Fig. ^h) with t = V = 1 and 
(a) A = 0, 5 and [7 = l7o = 0, (b) A = 5, [/ = 2, 5, 10, Uo = 0. 



discrete peak is slightly shifted through level repulsion. 

Results for similar parameters, but for the non-trivial 
topology of Fig. EJb), are shown in Fig. [5] Here the 
impurity forms a ring with its two neighbor O sites. The 
main difference is that the more complex "ring impurity" 
has more internal eigenstates, and these do not hybridize 
equally with the bath. For example, while there is a single 
eigenstates with one hole at the TM site, the isolated 3- 
site ring has 3 eigenstates with one hole on the ring. Two 
of these eigenstates involve the TM site (the third is the 
anti-symmetric combination of the O sites so it is a bath 
state). As a result, here we expect two one-hole continua 
in the AES two-hole spectrum, corresponding to states 
having one hole in one of these two ring eigenstates, and 
the second in the O bath. Similarly, while the TM atom 
has a single two-hole eigenstates, the 3-site ring has four 
such eigenstates in the Sz = sector. If Uo = 0, only 
3 of these involve the TM site, the fourth one is a bath 
state. As a result, for J7o = we should expect to see 3 
peaks marking this "multiplet" . 

A summary of the results is shown in Fig. |51 In panel 
(a), we plot Add{^^) for t = = 1, J7 = ?7o = and 
A = 0,5. The A = case (thin green line) looks very 
dilfcrcnt from its counterpart in Fig. [7{a): the two- hole 
continuum is located between [— 4f, 4t], as expected, but 
it has two broad resonances and there is a discrete peak 
just above it. This can be easily understood in terms of 
the ring "multiplet". For these parameters, the isolated 
ring supports three distinct two-hole eigenstates, with 
energies —2t,t,At. The two broad resonances are clearly 
associated with the former two. The third peak falls 
right at the two-hole continuum's upper edge, and level 
repulsion pushes it just above the band-edge. The one- 
hole continua overlap with the two-hole continuum, for 
these parameters, so they are not visible. 

For A = 5 (thick black line), the spectrum is rather 
similar to that shown in Fig. [TJb), despite the different 
topology. The reason is that for these parameters, the 3 
two-hole ring eigenstates are at — 10.63t, — 4t, 2.63t. The 
latter two fall inside the two-hole continuum, so only the 
former is expected to be visible. The result confirms this 
expectation. Similarly, the 2 one-hole ring eigenstates 
are at — 5.32t, 1.32t and should result in one-hole con- 
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FIG. 9. (color online). Impurity spectral function Ad^^dii^^) 
vs. CO for (a) the trivial topology and (b) the ring topology, 
with t = V = 1,A = U = 5 andUo =0,1,3. 



tinua spanning [— 7.32i, — 3.32<] and [— 0.68t, 3.32t], re- 
spectively. Again, only the former is (partially) distinct 
from the two-hole continuum, and indeed it is clearly 
visible in the spectrum. Thus, even though the multiplet 
structure is more complex for the ring topology, the re- 
sults do look quite similar to those for the trivial topology 
as most of the additional features fall inside the two-hole 
continuum. 

Similar conclusions are valid when we turn on the 
correlations, as shown in panel (b). The results are 
again qualitatively like those for the trivial topology, 
because the additional states in the ring multiplet fall 
inside the two-hole continuum and are washed out into 
broad, low-weight resonances. This is certainly true for 
U ^ t,\A\ ^ t, when the TM states become (to zero or- 
der in perturbation theory) eigenstates of the ring, with 
the remaining ring eigenstates involving only O sites and 
therefore being located inside the two-hole continuum. 
For finite values of U and A there are quantitative dif- 
ferences in the precise location of various features and in 
the shape of the continua. 

Finally, in Fig. [S] we briefly discuss the effect of the 
on-site repulsion at the O sites in the impurity approxi- 
mation. Panel (a) is for the trivial topology while panel 
(b) is for the ring topology. In both cases, t = V = 
1,U — A ~ 5, and three values Uq = 0,1,3 are shown. 
In the trivial topology, the main effect of Uq is to shift 
the one- hole continuum to lower energies by Uq, as ex- 
pected. This leads to a stronger hybridization with the 
discrete peak, which is consequently pushed to lower en- 
ergies as well. The weights in these features also very 
somewhat, except for the two- hole continuum which re- 
mains invisible on this scale. For these parameters, the 
results for both topologies are qualitatively similar for 
reasons explained above. 

To summarize, we find that the results for the two 
topologies are rather similar if |A| and/or U are large 
compared to t, V. If this is not the case their predictions 
can be quite different, however in such cases the impurity 
approximation itself becomes very suspect. 

When valid, the impurity approximation is fairly 
successful in reproducing the location of the multi- 
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FIG. 10. Momentum integrated AES spectral weight 
"^^k ^dd{k,ui) for the exact solution (full line) and our vari- 
ational approximation (dashed line) compared to the predic- 
tion of the impurity approximation with trivial topology (dot- 
dashed line), for i = 1/ = 1, [/ = A = 5, (7o = 0. 



plet peak, as shown in Fig. [TU] where we compare 
Addi^jj) for the trivial topology (dot-dashed line) against 
Add{k, uj) for the exact solution (full line) and our 
variational approximation (dashed line), for A — U — 
5, Uo = 0. The multiplet peak disperses very weakly for 
these parameters so all three curves show what looks like 
a discrete peak at low energies. The impurity approxi- 
mation is less successful in predicting the correct band- 
width for the continuum visible at higher energies; when 
analysing experimental data, this may result in wrong 
values assigned to the hopping parameters. In any event. 
Fig. [TO] confirms again that our variational approxima- 
tion is more accurate than the impurity approximation 
even for this simplest possible problem. 



V. DISCUSSION & CONCLUSIONS 

In this work we introduced a novel variational ap- 
proach to calculate two-particle Green's functions needed 
for AES spectroscopy. Here we applied it to a simple ID 
model with otherwise full bands, so that we can com- 
pare its predictions against the exact solution available 
for this case. We also calculated the results in the impu- 
rity approximation, for two different ways of connecting 
the impurity to the O bath. 

Both approximations are reasonable to use only if the 
single-electron parameters are such that one-hole eiges- 
tates with the hole primarily at the TM site (i.e., of TM 
character) are energetically well separated from those 
with the hole primarily in the O band. If this is not 
true, then the basis states projected out in these approx- 
imations are energetically close to those kept within the 
calculation, and the results are not sensible. 

If this condition is satisfied, we believe that these re- 
sults clearly show that our variational approximation is 
superior to the impurity approximation. Not only does it 
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produce momentum-resolved results, but the location of 
the continua that appear in the AES spectrum is in fair 
agreement with that expected from the convolution of the 
one-hole spectrum, unlike for the impurity approxima- 
tion. Since the one-hole spectrum can be obtained from 
ARPES measurements (and, from a theoretical perspec- 
tive, can be calculated with methods analogous to those 
we use here to find the two- hole propagator), a combina- 
tion of the two spectroscopies can be used to determine 
which features in the two-hole spectrum come from this 
convolution. Any other features must be due to on-site 
correlations, and can therefore help pinpoint the values 
of various on-site interaction energies. 

The method of Ref. [25| can be extended straightfor- 
wardly to compute Green's functions for problems with 
three or more holes, which would be relevant for AES is 
systems with partially filled bands, because the resulting 
equations of motion can still be grouped into simple re- 
currence relations. In practice, a variational approxima- 
tions along the lines we proposed here will significantly 
decrease the computational cost while removing from the 
spectrum only low-weight features whose location can be 
inferred apriorily from knowledge of the one-hole spec- 
tra, as already discussed. This opens an efficient way to 
make progress on understanding AES spectra for systems 
where the d orbitals of the TM are only partially filled, 
so long as the states of TM character are not too close 
to the O character states. Such work is now in progress. 
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Appendix A: Continued-fraction solution for the 
impurity approximation 

Since we lose the translation invariance, here we have 
to work explicitly in the real space. In the following we 
generalize the continued-fraction approach to the Ander- 
son impurity problem. 

For technical reasons which will become apparent soon, 
it is more convenient to group the TM together with the 
two O sites it directly hybridizes with into "site" 0, and 
to index the O atoms to its right/left as ±1,±2, etc. 
The two O inside "site" will be called "a" and "b", 
respectively. This indexing is shown in Fig. 1111 Amongst 
other things, this indexing makes it easy to study the 
two situations depicted in Figs. [TTJa) and (b), where the 
hopping between the two central O is turned off/on. 

We use the operator identity (a; + E^^ + z?; — 
'Hi)Gi{uj) ~ I to generate EOMs for the Green's func- 
tions. For example, for any \i\ > 1, |j| > 1, this 
results in: (w -|- + J7o(2 — Sij) + i'r])Gi^ji{uj) = 
tGi^.j+il{uj) + tGi-i^jl{u;)+tGt^j-il{uj) + tGi+i^jl{uj), 
where Gia-ja'i^) — {dd\Gi{u!)\ia, ja') . In other words, it 
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FIG. 11. Actual indexing of the O sites used in this work. 
The three central atoms comprise "site" 0. 



links together propagators where the holes are at a dis- 
tance = n apart, to propagators where the holes are 
nil sites apart. This allows us to rewrite these equations 
of motion as simple recurrence relations by grouping to- 
gether in a vector Vn all the propagators where the holes 
are n sites apart. Note that this is true for states with a 
hole at "site" and one at any site \n\ > 1, which also 
enter into Vn and only link to propagators in Vn±i- 
To be more precise, let: 



/ 



V° 



Gia:i+n,a' i^^) 
Gi-l^a;i+n-l,a' (w) 
Gi-2,a;i+n-2,a' (w) 



V 



be a vector of infinite length which contains all the prop- 
agators with the holes n sites apart. Of course, there are 
three entries replacing the (0, n) entry, corresponding to 
sites a,b,d, and similarly for (— n, 0). Further, we note 
that for n > 1, the ordering of the spins is preserved 
under hopping, since the left-most hole cannot pass by 
the right-most hole with one hop. Thus, we can write a 
recurrence relation for these vectors, for any n > 1, as: 



In 
7n 



an 






an 



n-1 
n-l 











yn 



(Al) 



where 7„, a„ and /3„ arc very sparse matrices, whose el- 
ements can be read off from the equations of motion as 
in the periodic case. Further, as before, the tridiagonal 
form of the recurrence relation admits a continued frac- 
tion solution for these vectors. The rest of the procedure 
is the same as for the periodic case. 

Truncation schemes of continued fractions: 
(i) For this method to work we need to truncate the 
continued fraction at a large inter-hole separation M, 
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both for the periodic and impurity cases. In the periodic 
system, we start out by creating two holes on a TM site. 
The two holes can now dclocalizc in the system. However, 
the broadening 77 introduces a finite lifetime ('^ ^/v) so 
that < k,Q^dd\Gp{Lo)\k,n^aP > — > as n — > 00. Thus 
for a large enough n = Af , we can truncate the continued 
fraction for the V-vcctors by setting Vm^+i to zero. 

For the impurity problem, this is justified because 
Gia-i+n,a-' (w) is the Fouricr transform of the amplitude 
of probability that the two holes move from the impu- 
rity (TM) site to the sites i,i + n within a time r. The 
larger n is, the less likely this process becomes, hence the 
smaller these propagators must be. This is certainly true 
for the energy ranges spanned by eigenstates that favor 
having the holes on the impurity (TM) site, since then 
they are unlikely to wander very far away from it. How- 
ever, this is also true even if the holes preferred the O 
bath. This is because the broadening rj is equivalent to 
introducing a finite lifetime ^ I/t] for the holes, so they 
cannot move arbitrarily far in a finite time r. Of course, 
in this latter case the appropriate cutoff M increases as 
r] decreases. In practice, for both cases, we increase M 
until the results become insensitive to further changes. 

(ii) A further truncation is necessary for the impurity 
problem. This is because, with the ordering of the prop- 
agators used for the Vn vectors, we have the additional 
complication that all these vectors, and therefore all the 
sparse matrices a„,/3„,7„, are infinitely-dimensional. In 
order to calculate the continued fractions, we need to 
truncate the size of these vectors, as well. The reasons 
discussed above justify doing this if either |j| » 1 and/or 
|i -l- ^ 1. We use the following truncation procedure: 
for a fixed separation n between the holes, the maxi- 
mum distance either of the two holes is allowed to travel 
away from the TM is Rc + n, i.e., we truncate V^''^ at 

G_R,,cr:i?,,+„,o-'(cj) at thc tOp and G-R^-n+l;a;~B^, + l.a'{^^) 

at the bottom. Again, Rc is increased until results be- 
come independent of its value. As a final comment, we 
note that there are other ways of grouping the propaga- 
tors into vectors so that thc equations of motion still lead 
to a simple recursive relation. Different schemes have var- 
ious computational advantages and disadvantages, but 
they all convergence to the correct answer if the cutoffs 
are sufficiently large. This converged result is equvalent 
to the exact solution computed by the Cini-Sawatzky the- 



ory for the Anderson impurity problem. 



Appendix B: Choice of parameters 

For simplicity, we discuss the impurity approximation 
first, and the full periodic system second. 

For the impurity approximation, if we ignore the TM- 
bath hybridization, V ^ 0, then if both holes are at the 
TM site the energy of the state is E'q — 2 A — J7; if one hole 
is at the TM site and one is in the O bath, the minimum 
energy of such states is Eq — A — U — 2t — Uq. Finally, 
if both holes are in the O bath, the minimum energy of 
such states is E^^ — At — 2Uo- As a result, thc state with 
both holes at the TM atom is favorable energetically if 
A > 2t+Uo and 2A+U > 4t+2t/o. Thc second condition 
is automatically satisfied if the first one holds, since the 
TM on-site interaction is repulsive, U > 0. 

The condition A > 2t + Uo implies that AES should be 
more useful for Mott insulators than for charge transfer 
materials 1^ To see this, consider one-electron removal 
states from the GS. If the hole is removed from the TM 
site, the energy of the state is Ef^ — A — U, while if it is 
removed from an O site, the minimum energy is EQ — 2t — 
Uq. The latter is energetically more expensive than the 
former if A > 2t -|- Uq- This makes sense, because if the 
material was a charge transfer insulator the additional 
holes would prefer to stay in thc O bath and AES would 
be less sensitive to the TM atom specific properties. 

The analysis for the periodic system is very similar for 
V 0, because the energy differences between states 
having (i) both holes at the same TM site, (ii) one hole 
at a TM site and one in the O bath, and (iii) both holes 
in thc O bath, are precisely the same as for the impu- 
rity case. Because of this and also in order to be able 
to meaningfully compare with impurity approximation 
results, we use the same parameters in both cases. 

Note that in thc periodic system we can also have (iv) 
the two holes at different TM sites. For V = 0, these 
states has energy Eq — 2A — 2U and are always energeti- 
cally favored compared to having both holes at the same 
TM site. Of course, these are the states projected out in 
our variational calculation. 
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